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Abstract. Halo mergers and shock waves play a crucial role in the process of hierarchical 
clustering. Hydrodynamical simulations are the principal investigation tool in this field for 
theoreticians, and predict that a by-product of cluster formation and virialisation is the injec- 
tion of turbulence in the cosmic flow. Here I will summarise results from a series of recent 
works focused on the main stirring mechanisms acting on baryons: minor and major cluster 
mergers, and curved shocks. Unresolved turbulence has been treated with the implementa- 
tion of a subgrid scale model. Recent simulations show that the production of turbulence 
differs significantly for the warm-hot intergalactic medium (WHIM) and the intra-cluster 
medium (ICM), because of different stirring mechanisms acting in the two gas phases. 
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1. Introduction 

The driving of turbulence in the intergalac- 
tic medium is a natural consequence of the 
hierarchical growing of the large-scale struc- 
ture. Interestingly, turbulence is both linked to 
the thermal properties of the cosmic baryons 
(as a by-product of the virialisation mecha- 
nism) and to the non-thermal diffuse cluster 
emission, because of the role which it is ex- 
pecte d to play in the acceleration of cosmic 
ravs (iFerrari et al.l l2008t iBrunetti & LazarianI 
l201lh and in the amplification of magnetic 
fields (ISubramanian et al.L l2006t iRvu et al.L 
120081) . 

The study of turbulence in the frame- 
work of the physics of galaxy clusters turns 
out to be challenging. From an observa- 
tional viewpoint, important progress has been 



made by measuring resonant scattering sup- 
pression (IChurazov et al.l 120041: IWerner et al.l 
2QQ2) and with XMM-Newton observations of 
cluste rs with a compact core (ISanders et al.l 
Numerical simulations, on the other 
hand, have to cope with the big scale separa- 
tion between the integral length scale for tur- 
bulence injection and the Kolmogorov scale, 
where the kinetic energy is dissipated by 
viscous effects. In grid-based hydrodynami- 
cal codes the use of adaptive mesh refine- 
ment (AMR) can only partially cure this is- 
sue; many efforts have been put in designing 
refinement crit e ria sui t able for turbulent flows 
( Schmidt et al., '2009'; lap ichino & NiemeveA 
2008; Vazzaetal., 2009). Alfliough recent 
simulations fin ally reach a fairly large dynami- 
cal range (e.g.. lVazza etaill2010l) . It IS notpos- 
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sible, even with AMR, to resolve the whole tur- 
bulent cascade down to the dissipative length 
scale. For this reason, turbulence subgrid scale 
(SGS) mo dels ('Sc annapieco & Bruggenil2008t 
iMaier et al.. 2009) provide the most physically 
motivated way of studying the effect of un- 
resolved turbulence on the system (see also 
ISchmidt & Federrathll2010l for a thorough dis- 
cussion on this point). 

In the following I will briefly summarise 
some recent results on the study of the main 
driving mechanisms of turbulence in the cos- 
mic flow at cluster scales: mergers and curved 
shocks. The work is based on hydrodynamical 
simulations perforrn ed with the ENZO code 
dO'Shea et al.Ll2005h . In fliis code framework, 
a subgrid scale model was implemented and 
coupled both with the equations of fluid dy- 
namics a t resolved l ength scales, and with the 
AMR CMaier et al.L |2009|) . The resulting nu- 
merical scheme has been called FEARLESS 
(Fluid mEchanics with Adaptively Refined 
Large Eddy Simulations), and combines the 
adaptive refinement of regions where turbu- 
lent flows develop with a consistent mod- 
elling of the unresolved turbulence. This tool 
has proved to be extremely useful in the 
modelhng of turbulent clumped flows, like 
those in ga laxy clusters . The reader is referred 
elsewhere dMaier et al.L l2009t llapichino et al.L 
I2OI ih for numerical details and tests of the 
SGS model and of FEARLESS. Recent devel- 
opments and improvements in the modelling 
of astrophysical turbulence are described by 
ISchmidt & FederrathI (I2OIOI) . 



2. Simulations of cluster mergers 

Galaxy clusters evolve mainly by accretion of 
smaller clumps; in this process, turbulent ki- 
netic motions are driven in the ICM. In case 
of minor mergers, the shearing instability de- 
velops at the interface between the ICM and 
the subcluster gas, resulting in the injection 
of turbulence in the region past the subclus- 
ter motion. The problem has been studied by 
means of high-resolution, idealised hydrody- 
namical simulatio ns (e.g., Heinz et al] 120031 : 
iTakizawa 2005: lapichino et al. 2008h . and is 



relevant for the physics of merger cold fronts 
dMarke vitch & Vikhlinin, 2007). 

Minor mergers have been also explored 
in full cos mological simulatio ns of cluster 
evolution. In iMaier et akl (|2009|) we used the 
FEARLESS approach to study the role of 
merger-induced turbulence in the ICM. Since 
the flow in the ICM is subsonic, we found that 
gas kinetic energy (either resolved or SGS) is 
just a minor part of the cluster energy bud- 
get. The production of turbulence is closely 
correlated with mergers: this is clearly visible 
in Fig. [U where a projection of the SGS tur- 
bulence energy (computed using FEARLESS) 
is compared with a volume rendering of the 
baryon density, showing the cluster substruc- 
ture. For example, one can observe the small 
subcluster immediately on the left of the clus- 
ter core, which is moving downwards around 
the centre and has stirred the ICM in its turbu- 
lent wake, as indicated by the large value of the 
SGS turbulence energy in that region. The dis- 
sipation of unresolved turbulence results also 
in a larger core temperature and a higher core 
entropy. 

The importance of mergers and merger- 
induced turbulence for the cluster energy bud- 
get is apparent from the stu dy of major merge r 
simulations performed by iPaul et al"] (1201 ll) . 
From the study of a sample of mergers, it was 
found that the ratio of the turbulent to total 
pressure in the cluster core is larger than 10% 
for about 2 Gyr after a major merger. The scal- 
ing of the turbulence energy with the cluster 
mass (Fig. |2]) in the sample is consistent with 
M^^^, which is the same scaling expected for 
the thermal energy in the self-similar model. 
This result highlights again that virialisation 
and turbulence injection are two faces of a 
same physical process, the hierarchical forma- 
tion of cosmological structure. 

3. Turbulence production in the 
intergalactic medium 

There are several mechanisms that are able to 
stir the baryons and inject turbulence in the 
cosmic flow at cluster scales. In the previ- 
ous Section, cluster mergers were presented. If 
the effects of galaxy motions in the IGM and 
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Fig. 1. Visualisation of a region of 12.8 Mpc centred on a cluster (M = 5.95 x lO'"* Mq A"' at z = 0), at 
z = 0.209. In the left-hand panel a volume rendering of the baryon density is reported, whereas in the right - 
hand panel a projection of the SGS turbulence energy is shown. Simulation details in iMaier et al.l ilOO^) . 
This is a snapshot from a full movie of the simulation, available at the website http://www.ita.uni- 
-heidelberg . de/"luigi/movies . html. 
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Fig. 2. The turbulence energy is plotted against the 
merger ma ss Mmem; for the merging clusters of the 
sample of iPaul et alj ( l201lh . The solid line is the 
best fit to the data points, corresponding to a scal- 
ing law oc M'-*, whereas the dashed line (with 
a slope of 1 .66) is the best fit computed by excluding 
an outlier cluster (the point at log M = 14.0). 



AGN outflows are neglected, a remaining stir- 
ring mechanism is provided by the baroclinic 
vorticity generation (Kang et al., 2007). In fil- 
aments and cluster outskirts the unprocessed 
gas is accreted and shock-heated in the form- 
ing structures, and the injection of turbulence 



is a by-product of this gas accretion at curved 
shocks. 

Bot h stirring mechanisms were studied in 
detail in llapichino et al.l (1201 ll) . by means of a 
FEARLESS simulation in a cosmological box 
with a volume of (100 Mpc /z"')^. The anal- 
ysis was focused on baryons with tempera- 
ture larger than lO'' K, and a distinction was 
done between gas in collapsing structures with 
baryon overdensity 6 > 10^ and less dense ma- 
terial. For consistency with previous studies in 
this field, it was chosen to refer to these two 
baryon phases as to ICM and warm-hot inter- 
galactic medium (WHIM), respectively. 

In Fig. [3] the time evolution of the specific 
internal and SGS turbulent energies is shown 
for both baryon phases. The production of tur- 
bulence has clearly a different redshift depen- 
dence in the ICM and WHIM. We argue that 
this difference is due to the mechanisms of tur- 
bulence generation that are dominant in the two 
baryon phases. 

In the ICM, turbulence is produced mainly 
by merger events: for this gas phase, the broad 
peak of et at redshift between 1.0 and 0.65 
is consistent with the major merger phase of 
halos in the mass range 
10^'^ Mq (for instance, iGiocoh et all l2007h . 
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Fig. 3. Time evolution of the mass-weighted aver- 
ages of specific internal (solid lines) and SGS tur- 
bulent (dotted lines) energies, for the two baryon 
phases under investigation. The two lines in the up- 
per part of the plot refer to the ICM, and the other 
two to the WHIM. The lines are scaled according to 
the factors in the legends, in order to be accommo- 
dated in the same plot. 

For the WHIM gas, residing in cluster out- 
skirts, smaller clumps and filaments, the evolu- 
tion of turbulence is related to the amount of ki- 
netic energy processe d by the external shocks 
(ISkilImanetal.L l2008'). where turbulence is in- 
jected by the baroclinic mechanism. For fur- 
ther details, and for a discussion about the 
effects of dynamical pressure support on the 
gravitational contraction of the ga s, we refer 
the reader to llapichino et al.l (1201 ll) . 
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